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Abstract 

It is widely believed that minimal supersymmetric SU(5) GUTs have been ex- 
cluded by the SuperKamiokande bound for the proton decay rate. In the minimal 
model, however, the theoretical prediction assumes unification of Yukawa couplings, 
Yd = Y e , which is known to be badly violated. We analyze the implications of this 
fact for the proton decay rate. In a consistent SU(5) model with higher dimen- 
sional operators, where SU(5) relations among Yukawa couplings hold, the proton 
decay rate can be several orders of magnitude smaller than the present experimental 
bound. 

PACS numbers: ll.10.Hi, 12.10.Dm, 12.60.Jv, 13.30.-a, 14.20.Dh 



1 Introduction 



Supersymmetric Grand Unified Theories (SUSY GUTs) [T] provide a beautiful framework 
for theories beyond the standard model (SM) of particle physics. They combine several 
attractive ideas, namely supersymmetry and unification of matter and interactions. A 
crucial prediction of SUSY GUTs is the instability of the proton [2 , and the long-lasting 
search for proton decay has put a strong constraint on unified theories. 

The simplest models are based on the gauge group SU(5). The SM particles can be 
grouped into two multiplets per generation, no additional matter particles are needed. 
Hereby, the down quark and charged lepton Yukawa couplings are unified. The GUT 
scale is set by the unification of the gauge couplings around 2 x 10 16 GeV in the Minimal 
Supersymmetric Standard Model (MSSM). 

SU(5) based models have been studied in great detail. Recently the simplest ver- 
sion, minimal supersymmetric SU(5) lj, was claimed to be excluded due to the Su- 
per Kamiokande bound on proton decay jSlHj- The exclusion of the "prototype" GUT 
model is an important result and it is worth analyzing the underlying assumptions care- 
fully. 

One ingredient is the sfermion mixings [3] which are essentially unknown and which 
are neglected in refs. 011] . Taking these mixings into account one can suppress the proton 
decay rate below the experimental bound [qT|o] . Another important question concerns the 
failure of down quark and charged lepton Yukawa couplings to unify. To our knowledge, 
all previous analyses assumed exact unification at the GUT scale, Y d = Y e , and then used 
the down quark matrix to study proton decay. The decay width, however, is strongly 
dependent on flavour mixing and there is no reason not to take, for instance, the lepton 
matrix instead. 

The failure of Yukawa unification can be accounted for by adding operators induced 
by Planck scale effects Since the GUT scale is only about two orders below the 
Planck scale, differences between down quarks and charged leptons can be explained by 
such operators. In addition, they also affect the proton decay operators. 

In this paper, we start with minimal supersymmetric SU(5) and discuss the influence of 
flavour mixing on proton decay. After that, we will study the impact of higher dimensional 
operators on proton decay. In particular, we consider two simple models where the decay 
rate is well below the experimental limit. 

The outline of the paper is as follows: After briefly describing the supersymmetric 
SU(5) GUT model (Section^ and analyzing the dimension five operators (Section OJ), we 
discuss the results of the different scenarios in Section 0] Important and clarifying details 
are given in the Appendices. 
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2 Supersymmetric SU(5) GUTs 



We start this section by briefly describing the minimal supersymmetric SU(5) GUT model 
PP. It contains three generations of chiral matter multiplets, 

10, = (Q,u c ,e c ), , 

and a vector multiplet A(2A) which includes the twelve gauge bosons of the SM and twelve 
additional ones, the X and Y bosons. Because of their electric and colour charges, the 
latter mediate proton decay via d = 6 operators. At the GUT scale, SU(5) is broken to 
Gsm = SU(3) x SU(2) x U(1)y by an adjoint Higgs multiplet £(24). A pair of quintets, 
H(5) and H(5*), then breaks Gsm to SU(3) x U(l) cm at the electroweak scale. The 
superpotential is given by 

W = \m tr £ 2 + \a tr £ 3 + \H{b*) (£ + 3a) H(5) 

+ i Y? 10, 10, H(5) + V2Y? 10< 5* H(5*) . ^ 

The adjoint Higgs multiplet, 

acquires the vacuum expectation value (VEV) 

(£) = ( xdiag(2,2,2,-3,-3), 
so that the X and Y bosons become massive, 

M v = M x = M Y = 5V2g 5 a , (2) 

whereas the SM particles remain massless. Here g§ is the SU(5) gauge coupling. The 
components £s and £3 of £(24) both acquire the mass 

M s ee M 8 = M 3 = \m , 

while £(3,2) and £(3* ,2) form vector multiplets of mass My together with the gauge mul- 
tiplets. Finally, the mass of the singlet component £24 is \vn. 

The pair of quintets, H(5) and H(5*), contains the SM Higgs doublets, Hf and Hf, 
which break Gsm, and colour triplets, H c and H c , respectively. To have massless Higgs 
doublets Hf and Hf, while their colour-triplet partners (leptoquarks) are kept super- 
heavy, 

M Hc = Mjj c = 5Xa , (3) 
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the mass parameters of H(5) and H(5*) have to be fine-tuned O(-) ~ 10~ 13 . This is 
the so-called doublet-triplet-splitting problem. As we will see below, RGE analysis gives 
constraints on the masses of the new particles. 

Expressed in terms of SM superfields, the Yukawa interactions are 

+ Qi Qj H c + Y$ Qi Lj H c + Y£ «° e J H c + Y% u° H c , 

where 

Y u = Ygq = Y ue = Y\ , (5) 
Y d = Y e = Y ql = Y ud = Y 2 . (6) 

In particular the Yukawa couplings of down quarks and charged leptons are unified. While 
rrib = m T can be fulfilled at the GUT scale, it fails for the first and second generation. 
This problem can be solved by adding higher dimensional operators due to physics at the 
Planck scale so that [S] 

= §m tr S 2 + \a tr £ 3 + b¥£±- + c^- . (7) 

Now the masses of S 3 and S 8 are no longer identical, which will affect the constraints 
on the leptoquark mass. Including possible couplings up to order 1/M P1 , the Yukawa 
interactions read 

W Y = ~ e abcde (y? 10? 10? H e + 10? 10? ^Hf + ft 10? 10? H d -3- 
4 V A^Pl ^Pl 

ya y c 

'PI Mpi 

Then the Yukawa couplings are given by 



+ v/2 ( V« ftf „ lOf 5-„ + ft? ff„ 10f= 5> c + ft? ftf„ 10? jf- 5' c 



Mpi Mpi 

Here a/Mp\ ~ (9(10~ 2 ), and S and A denote the symmetric and antisymmetric parts of 
the matrices, respectively. Thus the three Yukawa matrices, which are related to masses 
and mixing angles at Mz by the RGEs, are determined by six matrices. 
From eqn. Q one reads off, 

y y = 5-^-/i 2 . (10) 
M P1 y J 
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Figure 1: Proton Decay via dimension five operators: They result from exchange of the lepto- 
quarks followed by gaugino or higgsino dressing. 



Hence the failure of Yukawa unification is naturally accounted for by the presence of 
h 2 - Note that we do not need to introduce any additional field at M G ut to obtain 
this relation; it just arises from corrections 0(a/Mp{). Therefore we call this model a 
consistent supersymmetric SU(5) GUT model. In the minimal model, Y qq = Y ue = Y u ; 
furthermore, one usually chooses Y qi = Y ud = Y d . Note, however, that the choices Y qi = 
Y u d = Y e or Y q i = Yd, Y ud = Y e would be equally justified. As we shall see, this ambiguity 
strongly affects the proton decay rate. 

Finally, in general right-handed neutrinos can be added as singlets in SU(5) models. 
With lj = is?, the Yukawa interactions read 

Wy = Y$ h 5* H(5) + M ij h lj , 

where M is the Majorana mass matrix with eigenvalues O(Mgut)- 



3 Analysis of dimension five operators 

The evolution of the proton decay rate based on dimension five operators involves a 
number of parameters and assumptions which have changed in analyses during the past 
years. In this section we therefore list the main ingredients of our quantitative analysis. 
Technical details are given the Appendices. 

Integrating out the leptoquarks in eqn. Pjl. two dimension five operators remain which 
lead to proton decay (fig. [IJ 7J , 



M, 



He 



5 Y% Y q k r (Qi Qj) (Q* L m ) + YZ Y u k d m (v° e j) (ug <f mJ 



(11) 



called the LLLL and RRRR operator, respectively. The scalars are transformed to 
their fermionic partners by exchange of a gauge or Higgs fermion. Neglecting external 
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momenta, the triangle diagram factor reads, up to a coefficient K, depending on the 
exchange particle, 

/ mf^¥ m^¥ = (4^ /(M; mi ' m2) ' (12) 

with 

/(M;to 1) to 2 ) = ^^3^ ln ]^ - ln ]^J > ( 13 ) 

where M and to^ denote the gaugino and sfermion masses, respectively. 

As a result of Bose statistics for superfields, the total anti-symmetry in the colour 
index requires that these operators are flavour non-diagonal [8J. The dominant decay 
mode is therefore p — > Kv. Since the dressing with gluinos and neutralinos is flavour 
diagonal, the chargino exchange diagrams are dominant [HIEI]- The wino exchange is 
related to the LLLL operator and the charged higgsino exchange to the RRRR operator, 
so that the coefficients of the triangle diagram factor are 

ICllll = 2g 2 , 1C RRRR = yy' . (14) 

Here y and y' denote the corresponding Yukawa couplings (cf. fig. |l(b)| ) and g is the 
gauge coupling. 

The leading process p — > K + u is used in the analyses of Goto and Nihei and 
Murayama and Pierce |4 : to exclude the minimal supersymmetric SU(5) model. 



Calculation of the leading process 



The Wilson coefficients C 5 l = YqqYqi an d C 5 r = Y ue Y ud are evaluated at the GUT scale. 
Then they have to be evolved down to the scale M SU sy, leading to a short- distance 
renormalization factor A s . The sparticles are integrated out, as described above, and 
the operators give rise to the effective four-fermion operators of dimension 6. Now the 
renormalization group procedure goes on to the scale of the proton mass m p ~ 1 GeV 
leading to a second, long-distance renormalization factor A[. The factors are discussed 
in Appendix |D] 

At 1 GeV, the link to the hadronic level is made using the chiral Lagrangian method 
|11| I12|. In ref. j^j, the Amplitude for p — > K + v is given as 

2m r , 



A{K + v) 



PC 



usdu i _ /~iusdv 



3m, 



-D + 



PC 



udsu i _ r^udsv 



m 



1+7^(3^+^; 

im b 



+ aC 



dsuu 
RL 



m. 



3m 



B 



(3F-D) 



(15) 



where 
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• a and (3 are the hadron matrix elements ^3] 

au L (k) = e aPl (o 
(3u L (k) = e a/3l (o 

from which all other elements can be calculated. In our analysis we need |T3] 
(K + \(us R )d L \p) 
(K + \(ud R )s L \p) 

and a (3 for R <-> L . ul(\s) denotes the left-handed component of the proton wave 
function, /„. = 131 MeV the pion decay constant. It is known that \a\ ~ \(3\ [T3] . 
and different calculations give 0.003 GeV 3 < (3 < 0.03 GeV 3 . The latest evaluation 
was done by the JLQCD collaboration at 2.3 GeV obtaining a = —0.015(1) GeV 3 
and (3 = 0.014(1) GeV 3 J3|. The systematic uncertainties are large, and since we 
want to study whether the experimental limit excludes minimal supersymmetric 
SU(5) or not, we use the smallest value (3 = 0.003 GeV 3 ; 

• tub — 1150 MeV is an average baryon mass according to contributions from dia- 
grams with virtual S and A [TTj ; 

• D = 0.81 and F = 0.44 are the symmetric and antisymmetric SU(3) reduced matrix 
elements for the axial- vector current |15j : 

• the coefficients Cll an d Crl are related to the LLLL and RRRR operators as 
discussed in Appendix O and given by eqn. (jG.l|) . Moreover, they include the 
renormalization factors A s and A\ as well as the coefficient /C (|14j) . the triangle 
diagram factor (fT3*j) and the suppressing mass of the leptoquarks, M# c . 

The first line of eqn. (|15|) is related to chargino exchange as shown in fig. El This formula 
is given in ref. ^H]- The authors of ref. [3] add the third term due to neutralino exchange 
of fig. |7(b)| Here we also include the corresponding diagrams of fig. |7(a)[ 
Finally, the decay width is given by [T2"] 



r u r) u l 



p(k) 



Ku^ul p(k) 



a 2m p 
fn 3m B 
a 



Comparing the LLLL and RRRR contribution 



The RRRR contribution has been ignored for a long time. However, as pointed out 
by Lucas and Raby [TH], this operator gives a significant contribution in SUSY SO(10) 
models. The reason is that the Wilson coefficients and hence the LLLL contribution are 
proportional to = |(tan/3 + ^-g), whereas the RRRR contribution is proportional 
to (tan (3 + ^^g) 2 - Therefore the latter is dominant for large tan (3, which is naturally the 
case for SO(10) models. Here tan/? defines the ratio of the vacuum expectation values 
of the Higgs doublets Hf and Hf. Since the RRRR contribution is proportional to the 
Yukawa couplings, it is dominated by the third generation. As long as the top mass was 
believed to be less than 100 GeV, it could be neglected in the analysis. Then the decay 
width is given by the LLLL contribution and can be suppressed sufficiently by adjusting 
the phase matrix given in eqn. (jA.lj) . 

In ref. [3], the RRRR contribution was studied in the minimal SU(5) model. It was 
found that the total width is even affected for low tan f3 because the phase dependence of 
p — > K + P^ andp — » K + u T now differs, so both channels cannot be reduced simultaneously. 

Supersymmetric particle spectrum 

Looking at the dressing diagram we notice that by taking the sfermions to be degenerate 
at a TeV, the triangle diagram factor (|13J) is given by 

/(M;m)= M (m 2 — M 2 — M 2 In — . (18) 

J V ; (M 2 — m 2 ) 2 V M 2 J m 2 K J 

Therefore the sfermions are usually assumed to have masses of 1 TeV. An exception 
is often made for top squarks. Since the off-diagonal entries of the mass matrix are 
proportional to nit, the mixing is expected to be large, with at least one eigenvalue 
much below 1 TeV. In analyses, one typically uses 400 GeV, 800 GeV, or 1 TeV for m t ~. 
For the other sfermions, the mixings are neglected. The proton decay rate is further 
suppressed by light gauginos and higgsinos. Note that the experimental limit for charginos 
is m^± > 67.7 GeV [Hj. 

Since proton decay is dangerously large, also the decoupling scenario [TH] has been 
studied, where the scalars of the first and second generation can be as heavy as 10 TeV [3]. 
Here, proton decay is clearly dominated by the third generation. 

As already mentioned above, one can constrain the leptoquark mass Mh c by examining 
the RGEs for the gauge couplings; the details are given in Appendix |H1 This analysis has 
been done in the minimal model for a long time already, first at one-loop level then at 
two-loop accuracy because of the large top Yukawa coupling. The most recent calculation 
leads to the constraint 

3.5 x 10 14 GeV < M Hc < 3.6 x 10 15 GeV (90% C.L.) , (19) 
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with Mh c well below the GUT scale. 

This constraint depends on the Higgs representations. Other Higgs representations 
can be chosen as well which yield a higher leptoquark mass (cf. PHI)- Moreover, we 
already pointed out that M3 = Mg no longer holds in the consistent model. Then Mu c 
changes by a factor of (M 3 /M 8 ) 5 / 2 and we easily estimate that M 3 = 2M 8 is enough to 
raise the limit to Mqut- in our calculation, we therefore choose Mh c = 2 x 10 16 GeV. 

Finally, one can also define a quantity Mqut for which one gets [I], 

1.7 x 10 16 GeV < M GUT = (M£ M s ) 1/3 < 2.0 x 10 16 GeV (90% C.L.) , (20) 
in good agreement with the region of the gauge coupling unification. 



4 Minimal versus consistent models 

In this section we want to discuss the decay rate both in the minimal and in the consistent 
SU(5) model. The diagrams are given in AppendixO Finally, we calculate the decay via 
dimension six operators. 



4.1 Minimal model 

As already discussed in Section|21 we can choose Y d or Y e for Y qi and Y ud to calculate 
the proton decay amplitude. Since the Yukawa couplings of down quarks and charged 
leptons do not unify, this ambiguity cannot be resolved in minimal SU(5). Despite this 
fact, however, in all previous analyses the relations Y q i = Y ud = Yd have been used. 

As discussed in Appendix |XJ two physical bases are used to calculate the decay ampli- 
tudes, with either a diagonal up quark matrix or a diagonal down quark matrix [3]. 
Assuming 

Yqq Y ue Y u , Yql Y u d Yd , (^1) 

the Wilson coefficients at the GUT scale can be written as 

cl = KKi = ( v u p) (Vckm v d ) , 

C 5R = Y ue Y ud = (Dii K* KM )(P* V CKM V d ) 



in the former and 



Cl L = Y d qq Y q \ = (V C T KM PV U V CKM )(V 



' qq ql \ " <jkm~ - v-^ivi/ v- u,/ i (23) 

in the latter case. Here T> u and T>d are the diagonalized Yukawa coupling matrices evalu- 
ated from Y u and Yd, respectively, V CKM is the CKM matrix and P is the additional phase 
matrix as given in eqn. ()A.1|) . 
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Minimal model, Y ql = Y ud = Y d 



a> 
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Figure 2: Decay rate T(p — > K + u) as function of tan/3 in the minimal model with 
Y q i = Y u d = Yd- The experimental limits are given by SuperKamiokande experiment |2()|I21|. 



We choose the parameters in eqn. ((T7j) as described in Section El and vary tan/3. Since 
the decay width is proportional to tan/3, low values are preferred to obtain a small decay 
rate. On the one hand, the top Yukawa coupling becomes non-perturbative for low tan/? 
since h t ~ -^s. Hence, we start at tan/5 ~ 2.5. Fig. |2] shows the results of the following 
three cases: (i) all sfermions have masses of 1 TeV; (ii) m 4 ~ is changed to 400 GeV; (hi) 
decoupling scenario, where the scalars of the first and second generation have masses 
of 10 TeV. The dash-dotted line represents the experimental limit r = 6.7 x 10 32 years 
as given by the SuperKamiokande experiment [^IIlEjj the dotted line is the new limit 
t = 1.9 x 10 33 years j2I]. As in refs. [HHE], the amplitude is always above the experimental 
limit. 



Next, we study the case 



Y = Y =Y 

1 qq 1 ue 1 u j 



Yql = Y, 



ud 



(24) 



in order to illustrate the strong dependence of the decay rate on flavour mixing and 
therefore on Yukawa unification. The Wilson coefficients now read 



and 



Ct L = {V u P){MV e ) , 
C% R = (V u M*)(P*MV e ) 

C d 5L = (M T PV u M)(V e ) , 
C* R = (P*M*V e )(M T V u ) 



(25) 



(26) 
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Minimal model, Y ql = Y ud = Y„ 




tan (3 



Figure 3: Decay rate T(p — > K + u) as a function of tan/3 with Y q i = Y u d = Y e . The mixing 
matrix Ai is taken arbitrary or Ai = 1. 



where Ai = UlU e replaces the CKM matrix V CKM . Note that the mixing matrix in Y u or 
Y d (cf. eqs. and (fOjl ) is still given by V^ K m- Since 7^ Y" e , the masses and mixing 

of quarks and leptons are different and Ai is undetermined. 

We first ignore mixing, i.e. Ai = 1, and calculate the decay rate; the results are 
shown in fig. El Without mixing, only scalars of the first and second generation take part 
so that the decay rate can be reduced significantly in the decoupling scenario where the 
triangle diagram factor (|13|) changes by almost two orders of magnitude. 

Now we take Ai totally arbitrarily and minimize the decay rate. As can be seen in 
fig. El it is possible to push the amplitude below the experimental limit even for smaller 
sfermion masses. In the case mi = 400 GeV, this is only possible for small values of tan f3. 
The fact that a sufficiently low decay rate can be found illustrates the dependence on 
flavour mixing and therefore the uncertainty due to the failure of Yukawa unification. 



4.2 Consistent model 

In this case the coefficients of the operators can be derived from the superpotential (|8|). 

V -Y 9 a f s 1 ° f s 
qq ~ 1 Mp/ 1 2M P / 2 ' (27) 
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and 



Y ql = Y 2 + - 3^-h 2 



Mpi Mpi 



(2? 

y d = Y 2 + 2— /i! + 2—h 2 . 
ud 2 Mpi Mpi 2 



Note that Y g / — Y U( 2 = Y e — Y^, which means that Y q \ and Y ud cannot be zero at the same 
time. 

It is instructive to express these Yukawa matrices in terms of the quark and charged 
lepton Yukawa couplings and the additional matrices / and h (cf. relations eqs. ()A.5[) ). 

Y A = Y A - - f A 

2 Mpi (29) 

Yqi = Y e + 5— — hi , 

Yud = Y d + 5-^-/ii . 

iWpi 

If one allows the (3,3)-component of f\ and fa to be 0{^ 1 ) ^> 1, proton decay via 
dimension five operators can be avoided, for instance, by satisfying 

4 5cx 5 a 

so that both C 5L = Y qq Y q i and C 5R = Y ue Y ud vanish. 

Even if we restrict ourselves to 'natural matrices', i.e. couplings 0(1), we can con- 
siderably reduce the decay amplitudes by a suitable choice of matrices. In the following, 
we will illustrate this with two simple examples where either the RRRR or the LLLL 
contribution vanishes at the GUT scale. 

The first model is given by 

Y qq = Y ue = diag(0, 0, y t ) , 

Y ud = diag(0, y s - y^ y b - y T ) , (31) 
Y q i = diag(y e - y d , 0, 0) , 

where yj are the Yukawa couplings of the fermions at Mqut- In this model the RRRR 
contribution vanishes completely because C l 5 ^ m = Y£> e Y^™ is zero whenever a particle of 
the first generation takes part. But according to figs. |6(d)| and |7(b)| at least one particle 
of the first generation is needed. Furthermore, only the decay channel p — > Kv e remains. 
As requested, all matrix elements are 0(1) or smaller. 
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Consistent model I 




tan (3 

Figure 4: Decay rate T(p — * K + i>) as function of tan/3 in the consistent model. 

After RGE evolution by means of eqs. (jB.3|) and (jB.4j) . the simple structure of Wilson 
coefficients changes slightly, but the RRRR contribution and the decay channel p — > Kv^ 
are still negligible whereas p — > Kv T becomes dominant. Fig. H] shows the results for 
different sfermion masses. The decay amplitude is always well below the experimental 
limit, in the case = 1 TeV even more than two orders of magnitude. 

Now we turn to the second model where the LLLL contribution vanishes at Mqut; 

Y qq = Yue = diag(0, 0, y t ) , 

Y ud = diag(y d - y e , y s - y b ) , (32) 
Y ql = diag(0, 0, y T ) . 

Now C l ^ m = Y*i Yqi" 1 is only different from zero for i = j = k = m = 3 and the decay 
has to be non-diagonal. Only the RRRR contribution with a low absolute value remains. 
After renormalization, the RRRR contribution is still dominated by third generation 
scalars so that decoupling of the first and second generation does not change the result. 
The LLLL operator contributes only via p — > Kv T . 

As shown in fig. the proton decay rate is even smaller in this model. Furthermore, 
due to the smaller (3,3)-component of h\ compared to the first model, it can easily be 
used for higher values of tan f3. 
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Figure 5: Decay rate T{p — » in the second model. 

4.3 Proton decay via dimension six operators 

For completeness we include proton decay via X and Y bosons 
decay mode is p — > e + 7r°. The decay width is given by 



. The dominant 



9u 

My 



A 



2\2 



1 + (1 + lKdl 2 ) 



(33) 



The enhancement factor A contains both a short-distance contribution A SD between the 
SUSY-breaking and GUT scales and a long-distance contribution A LD = A\ between 
1 GeV and the SUSY-breaking scale. A SB splits into three parts according to the three 
gauge couplings: 



A SD = 


~ai{M z )~ 


23 
30 f>! 


~a 2 (M z )~ 


3 
2b 2 


~a 3 (M z y 




"5 




"5 




a 5 



4 
3b 3 



2.37 , 



(34) 



where the first part is an approximate calculation 

With the Super-Kamiokande limit r(p — ► e + 7r°) = 5.3 x 10 33 years and using 
a = 0.015 GeV 3 of ref. ^3], the mass of the heavy gauge bosons has to satisfy the lower 
bound 



M v > 6.8 x 10 15 GeV , 



(35) 



roughly half of Mgut fl2"Uj) . Since r oc M v , the proton decay rate for My = Mqut is far 
below the detection limit which can be reached within the next years. 
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5 Conclusion 



We have recalculated the proton decay rate in supersymmetric SU(5) GUTs. In particular, 
we have emphasized the strong dependence of the decay amplitude for flavour mixing. 

Minimal SU(5) GUT is inconsistent since the predicted Yukawa unification, Yd = Y e , 
is badly violated. A consistent supersymmetric SU(5) model requires additional interac- 
tions which account for the difference of down quark and charged lepton masses. Such 
interactions are conveniently parameterized by higher dimensional operators. 

We have shown that such operators can reduce the proton decay rate by several 
orders of magnitude and make it consistent with the experimental upper bound. We are 
not aware of a mechanism which would naturally lead to the required relations among 
Yukawa couplings. But, on the other hand, proton decay also does not rule out consistent 
supersymmetric SU(5) models. 
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A The SU(5) Yukawa sector and specific bases 

Minimal model 

In the minimal theory, the SU(5) Yukawa sector of the superpotential reads 

W Y = \ Y{ 3 10, 10, H{5) + V2 Y 2 lJ 10, 5* H{5*) . 

From the superpotential one can immediately conclude that Y"i is a symmetric complex 
matrix. With Y u = Y\ and Y d = Y 2 , the Yukawa matrices have the form 

Y U = U U V U PUZ , 
Y d = U d V d U\ R . 

Here, P is an additional phase matrix with det P = 1 which is usually parametrized as 

P = e iLp diag(e^ 13 ,e^ 23 ,l) . (A.l) 

These phases cannot be absorbed by field redefinitions 23 . The CKM matrix is then 
defined as 

V CKM = UlU d . (A.2) 
The most general weak basis transformation which leaves the interactions invariant is: 

lOi-io;- =£4,10, , 

5* -> 5*' = V.jb) . 

Then the Yukawa matrices transform like 

Y u — > U T Y u IA , 
Y d ^U T Y d V . 

The superpotential of the SU(5) Yukawa interactions expressed in terms of SM super- 
fields is given by eqn. Pjl. Transforming the singlets fields by $ — > W$$, the superpo- 
tential transforms like 

W Y = Q T {U T Y u U W u ) u c H f + Q T {U T Y d V W d ) d c H f + e cT ( Wj U T Y e V) LH f 
+ | Q T (U T Y qq U) Q H c + Q T (U T Y q i V) L H c 
+ u cT ( W„ T U T Y ue U We) e c H c + u cT ( W U T U T Y ud V W d ) d c H c . 
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There are two possible physical bases now, namely diagonal up quark and diagonal 
down quark matrices, which can be realized by a suitable choice of all transformation 
matrices. With Y qq = Y ue = Y u and Y q i = Y u d = Y d , the Yukawa interactions read 

W Y = Q T V u u c H f + Q T (V CKM V d ) d c H f + e cT V e L H f 

+ \ Q T (V u P) Q H c + Q T (V CKM V d ) L H c (A.3) 
+ u cT (V u V* KM ) e c H c + u cT (P* V CKM V d ) d c H c 

in the first and 

W Y = Q T (Vj KM V u ) u c H f + Q T V d d c H f + e cT V e L H f 

+ \ Q T (VLu V u P V* KM ) Q H c + Q T V d L H c (A A) 

+ u cT (V u V* KM ) e c H c + u cT (P* V CKM V d ) d c H c 

in the second basis. The former is used in ref. ^Uj, the latter in ref. 

In principle, these formulae are only valid for unbroken supersymmetry where one 
can use the same transformations for the fermions and their supersymmetric partners. 
Broken supersymmetry gives small corrections to these transformations 0. 



Consistent model 

Expanding the superpotential by higher dimensional operators, the identities (0) and (JB|L 

Y u = Y qq = Y ue = Y\ , Y d = Yqi = Y ud = y 2 ) 

at Mqut n o longer hold. Instead, one can derive the following relations between the 
matrices: 







Y 

qq 


Y 


_ 5 a f A 

2M P / 2 ' 








Y u - 


Y 




5 a 
4Mpi 






Y u - 


Y 




5 a 
4Mpi 


Y d 


-Y e 


= Y ud - 


Y qi 










lY d + 


lY 


= Y 2 3 a 
Mpi 


hi , 




-Y e 


= Y ud - 


-Y d 


= 5 J hl . 
M P1 





(/ 2 5 + 3/ 2 A ) 



(A.5) 



The antisymmetric part of / 2 is determined by the difference between Y qq and Y ue , then 
only symmetric terms of f\ and / 2 remain. 
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B Renormalization group equations 

Yukawa couplings 

The one-loop renormalization group equations, in the MS scheme, can be written for 
general Yukawa matrices |26] 



16/T 



16tt 



16tt 



• dY u 
dt 
2 dY d 
~df 

2 dYe 
dt 



T u - G u (t) + l(bY u Yi + c Y d Y} 
T d -G d (t) + ^(bY d Y} + cY u Yl 
Ye (r e -G e {t) + hY}Y)j , 



Y„ 



Yi 



(B.l) 



where t = log n/M z and the traces T u ,T d , T e are given by 

T u = t r ( 3 Y u yj + 3 a Y d Yj + a Y^ Y e ) 
T d = T e = tr(3aY u Y^ + 3Y d Yj+ Y$Y e ) . 



(B.2) 



The constants a, b and c as well as the functions G u (t), G d (t) and G e (t), are summarized 
in the tabled) 





SM 




MSSM 




(a, b, c) 


(1> 1> -I) 




(0, 2, 1) 




G u (t) 




8g!(t) 


g 9l{t) + 3 g*(t) + 


f 


G d (t) 






^^W + 3^ 2 2 (t) + 




G e {t) 






§</?(*) + 3 




b x 


l n +T6 m 




2 n + 




h 






2n + |m — 6 




b 3 


In -11 




2n - 9 





Table 1: Coefficients to (|B.1|) and ()B.2|I . The running gauge coupling constant at 1-loop is given 
by gf (t) = 9i(0)/ ^1 — fff (0) • The integers n and m stand for number of generations and 
Higgs doublets, respectively. 
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The equations for the Wilson coefficients read [3] 



lfor^C*" = (-8 gl -6gl- § gf) + C^ kl (Y d Yl + Y U Y^ 



dt 



d 



+ C£? kl {Y}Y e ) j + (Y d Yj + Y u Yj) + C^ m (Y d Yj + Y U Y^ 



IQ^jClf = (Sgl -Agl) Cg? + Cg" (2 YjY u ) 



+ Clt 1 h Y}Y d Y + Cgg" (2 Y e Y}) k m + C£T (2 



(B.3) 



(B.4) 



Gauge couplings, Constraint on Mh c 

In minimal supersymmetric SU(5), the RGE at one-loop level are given by [10J: 

1 



2tt 



2 1\ M s 



3\ A A 2 A 

5 / Alz My 5 M# 



a 2 - 1 (M z ) = a^ 1 (A) + i- 



2 13\, M s ,„ A A nl A 
n - log -± + 2n - 5 log — - 6 log — + 2 log — 

3 6 / Mx Mz My M 3 



a 3 - 1 (M z ) = a 5 - 1 (A) + 



2tt 



■\n-2 ) log^ + (2n-9)log— - -41og— - + 31og-^- + log A 



M 2 



M, 



where Ms is the SUSY breaking scale. Using the combinations 



(-a, 1 + So^ 1 - 2^ 1 ) (M z ) = — 

Ztt 

(5«r 1 - 3a, 1 - 2a- 1 ) (M z ) = ±- 

Z7T 



M S 12 / Mh 
-2 log log , 

B M z 5 B \ M z \M 8 



Ma 



Slog — + 12 log M| 



M 



Ho 



(B.5) 
(B.6) 



one can derive constraints on the products M^ C (M 3 /M 8 ) 5//2 and M v \/ M 3 M 8 = Mq UT . 
At two-loop level, there are no simple analytic relations any more. 

Taking M§ = M3 = Ms as it was done for the constraint in eqs. (fl"9|) and (|20jl . it 
simply reads M Hc and Mq UT = My Ms- 



C Diagrams 

Fig. |H1 lists the diagrams for the decay p — > fC + z/ with chargino dressing; they are related 
to the first addend in (|15jl. Those with right-handed fermions incoming can be divided 
into two groups depending on the dressing before (fig. |6(c)| ) or after the decay operator 
(fig. |6(d)| ); we therefore call them RRLL and RRRR diagrams, respectively. The latter 
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d, 



u 



-Vkl 



x w- 



d, 



>---x- 



• u 



Y qq (l,l)Y ql (j,i) (j + 1) Y qq (i,j)Y ql (2,k) (i / 2) 

(a) LLLL triangle diagrams 



d, 




rib, 



d. 







^ — 




c > 











-Vkl 



Y qq (l,j)Y ql (l,i) (j ? 1) Y qq (j,2)Y ql (i,k) (i ? 2) 

(b) LLLL box diagrams 

di. 



d< 




a 



d' 



H r 



>---x- 



' u 



Y qq (j,2)Y ql (i,k) (i ^ 2) Y qq (i,j)Y ql (2,k) (i / 2) 

(c) RRLL diagrams 



2' 



d' 











c > 









i/ 



d' 



H, 



c 



-Vki 



Y ue (i,t) Y ud (j, i) (j / 1) y„ e (i,o r ud (i, i) (j / i) 

(d) RRRR diagrams 
Figure 6: Diagrams with chargino dressing 
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d 







— ^ — 




: > 





-Vki 



d c 







— ^ — 


hp > 


f s 











-Vki 



u. 



Y qq (j,l)Y ql (i,k) (i^j) 
(a) LLLL diagrams 



Y qq (j,l)Y ql (i,k) (i^j) 
(b) RRRR diagrams 



Figure 7: Diagrams with neutralino dressing 

case is the only one related to the RRRR operator and C 5 r because there are no right- 
handed neutrinos in the model. As discussed in Sectional the dimension five operators 
are flavour non-diagonal, hence several diagrams are suppressed. 

By interchanging down and strange quarks as incoming and outgoing particles, we get 
the diagrams due to the second addend. The diagrams of the last one cannot be realized 
by chargino exchange, so we look at those with neutralino exchange that are given in 

fig. El 

The coefficients Cll and Crl used in eqn. fjTKj) then read 



C 



LL/RL 



M, 



He 



C 5 l/5R A s Ai f(M; mi, m 2 ) 



(C.l) 



where C$l = Y qq Y q i is related to Cll and C^r = Y ue Y u d to Crl- They are evaluated at 
the GUT scale and A s gives the correction due to running from GUT to SUSY breaking 
scale where K, is determined. In practice, the Wilson coefficients C5 are renormalized 
by means of eqn. (jB.3|) and (jB.4|) and evaluated at SUSY breaking scale. Finally, A\ 
describes the renormalization of the d = 6 coefficients down to 1 GeV. 



D Renormalization factors 

The renormalization factors are crucial for analyzing the proton decay. Since there is 
some discrepancy in the literature, we want to discuss them here in detail. 

As already mentioned in Section El there are two ranges for the renormalization, 
namely the short-distance between M GU t and M SU sy and the long-distance between the 
latter and the proton mass at ~ 1 GeV leading to the factors A s and A\. The former 
is highly dependent on the top Yukawa coupling y t and can therefore not be calculated 
analytically. 

The renormalization group effects in SUSY GUTs have first been discussed in ref. |27j . 
At that time, not only the high top mass was unknown (m t = 20 GeV was assumed), 
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but since there were no data at Mz, the values at 1 GeV were taken to calculate the 
decay rate. Hence the renormalization factors As and A^ were defined, which include the 
running factor of the Yukawa couplings from low to high scale. In this work, we use the 
Yukawa couplings at Mz and Msusy and evaluate their values at Mgut- These are taken 
as input parameters for the calculation, so our factors A s and A\ differ from As and Al 
in ref. (231101 • For the long-distance part, this discrepancy was stressed in ref. [2~H] . 

Because of the high top mass, A s cannot be solved analytically [TU] and depends on 
the related particles. Hence the Wilson coefficients are evolved down to Msusy by using 
eqs. f)B.3|) and (|B.4|) . For simplicity, M SUSY is identified with the electroweak scale, so A\ 
describes pure QCD renormalization down to 1 GeV 1 , 



had/ 



a 3 (M z 



had J 



a 3 (m c 



_a 3 (m b )_ 



a 3 (m z ) 



2:i 



1.43 



(D.l) 



with /i had = 1 GeV. 
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